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Abstract. Spectroscopic observations of 48 emission-line objects of M 33 have 
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been obtained with the multi-object, wide field, fibre spectrograph AF2/WYFFOS 
at the 4.2m WHT telescope (La Palma, Spain). Line intensities and logarithmic 
extinction, c/j, are presented for 42 objects. Their location in the Sabbadin & 
D'Odorico diagnostic diagram (Ha/[S ill vs HWfN ill) suggests that >70% of 
the candidates are Planetary Nebulae (PNe). Chemical abundances and nebular 
physical parameters have been derived for the three of the six PNe where the 
4363A [O in] emission line was measurable. These are disc PNe, located within 
a galactocentric distance of 4.1 kpc, and, to date, they are the farthest PNe 
with a direct chemical abundance determination. No discrepancy in the Helium, 
Oxygen and Argon abundances has been found in comparison with corresponding 
abundances of PNe in our Galaxy. Only a lower limit to the sulphur abundance has 
been obtained since we could not detect any [S ill] line. N/H appears to be lower 
than the Galactic value; some possible explanations for this under-abundance are 
discussed. 
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1. Introduction 

Chemistry of galaxies has been studied for a long time using spectroscopy of the inte- 
grated stellar light. There are, however, limitations to this method, such as the large 
abundance variations from star to star and the consequent difficulty in the interpreta- 
tion of the data (Mc William 1997), or the impossibility of using this technique in poorly 
populated regions, as in the haloes of elliptical galaxies. 

Planetary nebulae (PNe) with their prominent emission lines can be used to measure 
the metallicity in nearby galaxies. In spite of the importance of the topic, few spectro- 
scopic data have been obtained for PNe in external galaxies, except for the Magellanic 
Clouds. PNe in the LMC and the SMC have been extensively studied by several authors, 
for instance Dopita et al. I|1997(l . Reyes et al. ©97), and Richer (1993|). To date, there 
are only few other galaxies. In fact, abundance analysis requires the determination of 
electron temperature and density, as well as the intensities of various lines of different 
ionization states. The most important line for temperature determination, [O in] 4363A 
is, for instance, from 50 to 200 times weaker than [O ill] 5007A (e.g. Ford et al. 2002 ). 
Before our work, the farthest galaxy with direct abundance measurements in PNe was 
M 31 (750 kpc, Freedman et al. I2001|l . Chemical abundances of PNe in M 31 have been 
investigated by Jacoby & Ciardullo (1999, hereafter JC99) and Stasinska et al. I|1998[) . as 
well as by Hyung et al. H2000|) who used the data from the previous work of Stasinska and 
collaborators. A total of 45 PNe have been studied. In the companion galaxies of M 31, 
spectrophotomctric data have been obtained for PNe of NGC 185, NGC 205 (Richer & 
McCall [THUS) and of M 32 (Stasinska et al. ITMSl Hyung et al. FMRty . In the Fornax 
galaxy, one PN is known and its chemical abundances have been derived by Danziger et 
al. (1978) and Walsh et al. (1997). Spectrophotometry of the two PNe in the Sagittarius 
dwarf elliptical galaxy was done by Walsh et al. I|1997(l . Chemical abundances of PNe 
in NGC 6822 were computed by Richer & McCall ifTMB) . Spectra of 5 PNe in the giant 
elliptical galaxy NGC 5128 (Centaurus A) at the distance of 3.5 Mpc were analyzed by 
Walsh et al. (1999). These are the farthest PNe where chemical abundances have been 
measured. This was done using diagnostic line ratios of various ions without, however, a 
direct determination of the electron temperature. 

The galaxy M 33 (NGC 598) is the third-brightest member of the Local Group. Its 
large angular size (optical size 53'x83', Holmberg 1958) and its intermediate inclination 
i=56° (Zaritsky ©99} make it particularly suitable for studies of spiral structure and 
stellar content (van den Bergh 2000). A total of 131 candidate PNe are known in M 33 
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distributed throughout the whole galaxy (Magrini et al. 120001 I2001afl . The objects were 
identified with the following two criteria: i) they should appear both in the [O in] and 
Ha+[N n] images but not in the continuum frame, and ii) they should have a stellar 
point spread function. In this paper, we present spectroscopy of 39 PN candidates and of 
9 unclassified emission- line objects discovered in our previous survey (Magrini et al. 120001 
hereafter MOO). In Section 2 we describe the observations, the reduction procedures and 
the flux measurement with their de-reddening. The analysis of the spectra and the deter- 
mination of the physical and chemical properties of the PNe are presented in Section 3. 
A critical discussion on the chemical abundances is given in Section 4. Conclusions follow 
in Section 5. 



2. Observations and data analysis 

Thirty-nine candidate PNe and nine unclassified objects with emission lines and a non- 
negligible continuum, selected from the list of MOO, were observed on October 16-17, 2001. 
We have used AF2/WYFFOS, the multi-object, wide field, fibre spectrograph working 
at the prime focus of the 4.2 m William Herschel Telescope (La Palma, Spain). The 
WYFFOS spectrograph was used with a single setup: the R600B grating (600 line mm -1 ) 
providing a dispersion of 3.0 A/pixel. The resulting spectral range, from 4300A to 7380 
A, included the basic lines needed for the classification of the objects as PNe and the 
determination of their chemical abundances. However, due to the arrangement of fibres 
at the spectrograph entrance, not all the spectra start at the same wavelength and for 
15 PNe the A4363A spectral region was not observed. The spectrograph WYFFOS was 
equipped with a 1024x1024 TEK CCD. We used the Small Fibre module which is made 
of 150 science fibres with 1.6 arcsec diameter (90 /xm) projected on the sky. 

Our targets, their astrometry and their absolute [O in] 5007A and Ha + [N n] fluxes 
came from our previous INT+WFC images (MOO, Magrini et al. 2001a). The accuracy 
in their positions is better than 0.5 arcsec r.m.s.. Our set of targets was chosen among 
the PNe spanning a range of 3 mag in [O ill] A5007 A from the bright cutoff of their 
luminosity function, and included also several of the unclassified objects of our sample, 
i.e emission-line objects with a non-negligible continuum emission (MOO). Targets are 
distributed rather uniformly over the face of the galaxy (Figure 1). We used the remaining 
fibres to take simultaneous sky/background spectra. A total of 13 science exposures of 
2400 s each were taken (5 in the first night and 8 in the second night) through light cirrus 
with 1.5 arcsec seeing. 

Several offset sky exposures using the same fibre configuration were taken before and 
after the M 33 observations in order to do a correct sky subtraction. 
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Fig. 1. Target positions in the M 33 galaxy. The i.d. numbers are from MOO. The emission- 
line + continuum objects are shown with their i.d. number, preceded by a capital C. The 
image is from the ESO-Digitized Sky Survey and its size is 40' x 40 '. North is at the 
top, East to the left. 

2.1. Data Reduction 

The data were reduced using the IRAF multi-fibre spectra reduction package DOFIBER. 
The sky subtraction represented the most difficult part of the reduction. As already 
mentioned, we used some of the fibres to monitor the sky background during the science 
exposures and we also had several offset sky exposures with the same fibre configuration. 
We used the sky fibres in our exposures of M 33 to monitor the relative intensity of 
the atmospheric emission lines, which are known to vary during the night. However, 
these sky fibres lie in different positions in the spectrograph focal plane and consequently 
produce spectra with different spatial and spectral instrumental profiles. For this reason, 
they cannot be used for a direct sky subtraction. There are instead no variations in the 
instrumental profiles for each individual fibre between the offset sky frame spectra and 
the science frames, as fibres are fixed at the spectrograph entrance. We therefore used 
these offset sky spectra to do the sky subtraction, after correcting the relative intensity 
of the atmospheric lines using the mean sky spectrum computed from the sky fibres 
observed at the same time as the science targets. 

2.2. Flux calibration 

Relative data flux calibration was obtained taking spectra with several fibres of the 
spectrophotometric standard stars G191B2B fOke ll990[> and 40 Eri B fOke ll974[> . Using 
AF2, the absolute flux calibration cannot be done, since one is limited by the fixed fibre 
diameter. In addition, it is impracticable to observe a spectrophotometric standard star 
in all the science fibres in order to have a sensitivity function for each fibre. In order 
to test whether a mean sensitivity function can be used for all the fibres, the standard 
star G191B2B was observed in nine fibres and the sensitivity functions of those fibres 
were analyzed. The variations of the spectral response of each fibre are within ^1% in the 
whole spectral range, as shown in Figure 2. This result allows us to use a mean sensitivity 
function for all the fibres. 

2.3. Flux measurement and de-reddening 

Emission line fluxes of the spectra were measured using the MIDAS/ ALICE package. We 
also used the IRAF task TWOFITLINES (IACTASKS, Acosta HMQll to estimate errors 
on the fluxes. This task computes the uncertainty of model parameters by a number 
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Fig. 2. The percentage variation from fibre to fibre of the sensitivity function. 
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of repeated fits to simulated data, generated from the best model with random noise 
added. The absolute errors on the measurement of the emission-line strengths are given; 
percentage errors (not including flux calibration errors) on the H/3 line are shown in 
Tabled The background noise and the sky subtraction are the greatest sources of error 
and they affect uniformly the whole spectral range where emission lines were measured. 
Thus, the absolute error of each emission line of our spectra is approximately constant 
and equal to the absolute error of H/3 line. Moreover, systematic errors due to the flux 
calibration, including the errors on the Oke fluxes (~ 0.01 mag), were estimated to be 
approximately 5-10% along the whole spectral range. 

The observed line fluxes have been corrected for the effect of the interstellar extinction. 
The extinction law of Mathis (|1990fl with Ry=3.l has been used. The eg, which is the 
logarithmic difference between the observed and un-reddened H/3 fluxes, was determined 
comparing the observed Balmer I(Ha) /I(H/3) ratio with its theoretical value. All fluxes 
are presented on a scale relative to H/3, where F(H/?) is set to 100. Observed fluxes and 
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cp for the objects in whose spectra we could extract information are shown in Tabled 
The percentage errors on the H/3 fluxes are presented in the last line of Table ^ 

3. Analysis of the spectra 

Out of the 48 targets observed, the spectra of 3 candidate PNe and 2 unclassified ob- 
jects were so weak that no usable information could be extracted, whereas one of the 
unclassified objects, namely MCMPc9, turned out to be a normal star. The nature and 
properties of the remaining 42 objects are discussed in the following. 

3.1. The Balmer decrement: PNe extinction 

The weighted mean value of the logarithmic extinction, c@, of the PN candidates (Table 1) 
is 0.73±0.21. Approximately 30% of the candidates have a very small cp and presumably 
they lie in front of the galaxy. The other PNe have cp greater than 0.5: they lie on the 
far side of M 33 and/or in regions of the spiral arms where the extinction is higher. The 
mean logarithmic extinction of our disc PNe is similar to that found by JC99 in three 
PNe on the disc of M 31 (cp= 0.50; JC99). 

3.2. Diagnostic diagrams 

Magrini et al. (2000) introduced the excitation parameter R= gj^njj^ as a crude indi- 
cator of the nature of extragalactic emission line objects found in narrow-band surveys, 
arguing that most objects with R > 1 (and point-like at the distance of M31 and M33) 
are expected to be PNe. Ciardullo et al. (2002) suggest a slightly higher limit (R > 1.6) 
for PN candidates in the brightest 1-magnitude bin of the PNLF. Nevertheless, using 
R > 1.6, we would exclude a large number of low excitation PNe. In fact about 40% 
of confirmed Galactic PNe have R < 1.6 (see MOO, Figure 3). In the following analysis 
on the nature of candidate PNe, we consider adequate for PNe values of R from 0.3 to 
larger than 1.6, as expected for a 'normal' (i.e. Galactic) population of PNe. However, 
to identify bona fide PNe, more precise indicators should be used when, as in our case, a 
larger part of the spectrum is available. The empirical emission line diagrams (Sabbadin 
& D , Odorico ll976l Garcfa-Lario et al. 119911 Riesgo-Tirado & Lopez 2002J have proved to 
be very useful in distinguishing among different classes of Galactic ionized nebulae, the 
most effective one being probably that displaying the Ha/[N n] vs. Ha/[S n] line ratios, 
as it neatly separates PNe from H n regions and supernova remnants. This diagnostic 
diagram is especially valuable for extragalactic objects since it involves relatively strong 
lines and avoids information less easily obtainable in the UV or IR spectral ranges. We 
have recalculated the PNe zone in this diagram using the most precise Galactic and extra- 
galactic data available: Kingsburgh & Barlow 1)19941 hereafter KB94), Kaler et al. {T996), 
Cuisinier et al. (2000) for Galactic PNe; de Freitas Pacheco et al. I|1993[) . Vassiliadis et 
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Identification 




MCMP8 


MCMP9 


MCMP17 


MCMP18 


MCMP23 


MCMP25 


MCMP28 






PN 


PN 


- 


PN 


H II r 


PN 


PN 


4685.7 


He II 


- 


- 


- 


- 


- 


- 


29.1 


4861.3 


H/3 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


4958.9 


[O in] 


160. 


467. 


206. 


255. 


67.3 


221. 


396. 


5006.8 


[O in] 


454. 


1275. 


469. 


641. 


165. 


803. 


1198. 


5875.7 


He 1 


- 


33.9 


- 


- 


22.0 


- 


13.6 


6562.8 


Ha 


636. 


741. 


285. 


350. 


364. 


446. 


509. 


6583.4 


[N 11] 


52.3 


108. 


- 


75.3 


29.8 


61.0 


45.3 


6678.1 


He 1 


14.8 


20.8 


- 


20.3 


- 


- 


11.5 


6716.5 


[S ii] 


23.9 


17.5 


18.3 


65.0 


70.6 


132. 


50.2 


6730.8 


[S ii] 


23.9 


21.6 


18.8 


39.2 


50. 


- 


36.6 


7065.3 


He 1 


23.1 


- 


- 


35.0 


- 


- 


- 


7135.8 


Ar III 


62.7 


80.2 


- 


62.8 


- 


- 


48.9 


7325 (*) 


[O II] 


22.7 


15.2 


- 


- 


- 


- 


19.0 


c fi 




1.15 


1.36 


0.00 


0.30 


0.35 


0.64 


0.83 


H/3 error 




2% 


10% 


15% 


4% 


8% 


50% 


6% 






MCMP31 


MCMP38 


MCMP40 


MCMP41 


MCMP42 


MCMP45 


MCMP49 






H II r/PN 


PN 


PN 


SNR 


SNR 


PN 


PN 


4685.7 


He II 


13.6 


- 


- 


30.6 


- 


- 


16.5 


4861.3 


H/3 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


4958.9 


[O in] 


197. 


134. 


517. 


119. 


339. 


218. 


208. 


5006.8 


[O in] 


314. 


158. 


980. 


265. 


840. 


423. 


487. 


5875.7 


He I 


15.4 


91.9 


- 


64.6 


- 


14.6 


21.4 


6562.8 


Ha 


288. 


285. 


662. 


736. 


329. 


296. 


595. 


6583.4 


[N II] 


103. 


71.9 


964. 


250. 


130. 


96.1 


88.8 


6678.1 


He I 


13.3 


- 


- 


- 


- 


- 




6716.5 


[S II] 


41.0 


28.4 


183. 


189. 


284. 


26.6 


53.0 


6730.8 


[S II] 


30.7 


- 


135. 


158. 


123. 


23.6 


32.7 


7135.8 


Ar III 


32.6 


- 


- 


75.1 


- 


15.9 


- 


cp 




0.01 


0.00 


1.20 


1.36 


0.20 


0.05 


1.05 


H/3 error 




8% 


15% 


40% 


20% 


50% 


10% 


8% 






MCMP53 


MCMP60 


MCMP61 


MCMP65 


MCMP68 


MCMP69 


MCMP71 






SNR 


PN 


H II r/PN 


PN 


PN 


PN 


PN 


4363.2 


[O in] 


- 


16.4 


- 


20.2 


- 


- 


- 


4685.7 


He II 


- 


- 


78.2 


- 


- 


- 


- 


4740.2 


Ar IV 


- 


- 


- 


- 


- 


13.4 


- 


4861.3 


H/3 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


4958.9 


[O in] 


- 


568. 


274. 


510. 


160. 


288. 


365. 


5006.8 


[O in] 


- 


1475. 


637. 


1259. 


520. 


748. 


553. 


5875.7 


He I 


- 


- 


- 


- 


32.4 


41.9 


- 


6562.8 


Ha 


1177. 


484. 


320. 


567. 


365. 


555. 


658. 


6583.4 


[N II] 


428. 


221. 


98.1 


163. 


134. 


693. 


539. 


6716.5 


[S II] 


410. 


54.1 


43.3 


71.7 


67.2 


201. 


53.7 


6730.8 


[S II] 


356. 


48.6 


49.7 


74.8 


63.3 


163. 


63.0 


7065.3 


He I 








22.5 


16.2 


68.2 




7135.8 


Ar III 




52.0 


84.8 


119. 


39.1 


145. 




7325(*) 


[O II] 




29.6 




40. 








cp 




2.03 


0.76 


0.17 


0.97 


0.35 


0.95 


1.19 


H/3 error 




50% 


20% 


20% 


10% 


8% 


20% 


50% 



Note: (*) Sum of the [O n] doublet lines at 7319.6 and 7330.2 A. 
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Table 1 - continued 
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±vx ivi r i j 


MPMP77 
±vx v_,' ivi r i I 


ivi ±vx r / o 


MPMPQ1 
ivi ivx jt y x 


MPMPQ3 

IVI IVX JT y o 


MPMPQfi 

IVX IVI JT c/ U 


MPMP1 01 

1V1\_^1VXJ^ xux 






PN 




H II r 


PN 


PN 


PN 


PN 


A 1(^1 9 


IP> tttI 
[U IIIJ 








Q fi 

y .o 


fi nn 
o.uy 




inn 

1U.U 


A A 71 K 
4:4 I 1.0 


He n 






7 ^ 
( .00 


12 3 


1/17 
14. ( 






4685.7 


XX „ j j 

llC 11 








ICQ 

lo.o 


O.OZ 






4740 2 




11 9 














A Sfil "3 
4001.0 


XX Q 

rip 


i nn 

1UU. 


i nn 

1UU. 


i nn 

1UU. 


i nn 

1UU. 


i nn 

1UU. 


i nn 

1UU. 


1 nn 

1UU. 


4yoo.y 


ro tttI 

[U IIIJ 


A ^ 1 
401. 


ZOO. 


51 1 


448 


4:U0. 


404 


40U. 


DUUO.O 


ro tttI 

[U IIIJ 


1212 


f oo. 


1 ^7 
10 I . 


1 1A a 
low. 


i 9nn 

1ZUU. 


1 1 QQ 

i iyy . 


1 Qar; 
I00O. 


Oo /O.i 


He I 






OQ Q 

zo.y 






1 Q Q 

iy .y 


24 3 


6562 8 


Hct 


366 


600 


427 


433 


430 


438 


499 


6583 4 


TN tt! 
[IN 11J 


151 




73 2 


153 


91 8 


96 6 


74 1 


6678 1 


He I 


17 3 




8 66 


11 


5 95 


10 1 


11 4 


fi71 fi ^ 
1 1 o . 


\Q. tt! 
[d IIJ 


10 4 




^n i 

OU. 1 


Q A 

y .4 


1Q C 

iy .o 


n 

OO.U 


15 2 


6730 8 


ttI 


18 7 




35 2 


13 2 


21 6 


22 5 


11 3 


7065 3 


He I 


24 6 




7 32 


10 5 




23 




7135 8 




45.3 




23 7 


31 2 


27 3 


36.5 


38.6 


7319.6 


ro iii 






12.1 


10.5 


10.8 




12.3 


7330.2 


ro iii 






10.9 


7.7 


8.3 




10.8 


C(3 




0.36 


1.06 


0.58 


0.60 


0.59 


0.62 


0.80 


H/3 error 




8% 


50% 


9% 


2% 


5% 


5% 


3% 






MP MP 1 0^ 

IV X IVX l l-\JO 


MPMP1 OR 

IVI ±VX I 1UU 


MPMP1 99 

IVI IVX X^ X ^ 


MPMP1 9^ 

IVX IVI JT X O 


MPMP1 97 

IVX IVI X li ( 


MPMP1 "3D 

IVIV^IVXX^ 1(>U 


MPMP1 ^9 

IVI IVX I IOa 






PN 


PN 


H TT r/PN 

XX 11 1 / JT IN 


PN 


H II r 




H TT r/PN 

XX 11 1 / X IN 




ro tttI 

[V IIIJ 






1 1 fi 
l.OO 










A ATI ^ 

4:41 f 1.0 


He II 






S fiS 
o.Oo 






9n k 

ZU.O 


7 QQ 

1 .yo 


/I Sfil 1 


lip 


1UU. 


i nn 

1UU. 


100 


100 


i nn 

1UU. 


i nn 
1 U u . 


100 


4958 9 


ro tttI 

[U 111J 


367 


228 


58 


189 


61 8 


50 9 


90 8 


OUUD.O 


rn tttI 

[U 111J 


787 


265 


179 


583 


194 


124 


234 


Oo 1 O. / 


He I 






25 4 


14 9 


97 n 
Z ( .u 


ion 
lo.U 


ZO.O 


fi^fi9 S 
DOOa.o 


na 


OOZ. 


( OO. 


fiSQ 

ooy . 


4:00. 


7^ ^ 
/ OO. 


oyo. 


621 


00oO.4 


rw ttI 

[IN IIJ 


A 1 * fi 


A 7fi 
4 / O. 


101. 


123 


221 


i fin 

10U. 


ss 7 
00. / 


6678 1 


He I 






10 9 




77 3 




21 


fi71 fi ^ 
O i 1 . 


rc; ttI 
[q IIJ 


18 2 


SI "3 
ol.O. 


fil fi 
Ol.O 


Oo.o 


148 




/in s 

4U.o 


O ( oU.o 


rc; ttI 
[fc> IIJ 




^ 1 
OO. 1 


OO.o 


oo.y 


1 1 Q 

i iy . 




Tfi Q 

00. y 


7065 3 


He I 






10 5 








21 7 


7135 8 




15 3 




29 7 


34.1 


95 5 




50.4 


7319.6 


ro iii 






10.0 










7330.2 


[O III 
[v nj 






10.3 


26.9* 








c f3 




0.42 


1.45 


1.26 


0.68 


1.39 


0.47 


1.11 


H/3 error 




12% 


40% 


2% 




12% 


10% 


4% 






MPMP1 "34 

IVX V_> IVX X^ li)^ 


MPMPr^ 

IVX IVI JT X. U 


IVX v_.' IVI JT X \_) 


MPMPrl S 
ivi ivx r lio 


MPMPrl 

IVX v_.' IVI X LOO 


MPMPr4S 

IVX v_. IVI JT X 4: 


MPMPr^D 
ivi ivx r lju 






PN 


H II r 
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[S II] 


51.2 


56.4 


29.2 


99.0 


29.0 
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Ho I 
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1.36 


H/3 error 
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Note: (*) Sum of the [O 11] doublet lines at 7319.6 and 7330.2 A. 
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Fig. 3. Diagnostic diagram (originally from Sabbadin & D'Odorico. ll976fl . The Galactic 
H II regions and SNRs areas are from Garcfa-Lario et al. ITMT1 The PNe region is 
drawn using recent data of Galactic and extragalactic PNe (see text). Our candidate 
PNe are represented by triangles (R>1.6), circles (0.3<R<1.6) and squares (R<0.3). 
The emission-line objects with non-negligible continuum are represented by stars. 




log(Halpha/[SIlJ) 

al. lfT5§2)> and Stanghellini et al. fZHtf^ for LMC and SMC; JC99 and Stasinska et al. 
(1998 ) for M 31; Walsh et al. I|1997jl for Fornax. Note that the new limits shift the upper 
border of the PNe band more into the H n region box. 

The M 33 PN candidates were put in this diagram (Figure 3) coded by their R 
parameter. In this way we supplement the diagram with the information of the [O in] 
line which is the most important optical line characterizing a PN. We analyze in the 
following the location of all objects in this diagram. Five of the nine objects with non- 
negligible continuum (see MOO) have good data for the involved lines (Table 1). Two 
of them (MCMPcl8 and 43) are near the SNRs locus; both have strong [S n] lines (in 
MCMPcl8 the [S n] lines are stronger than [N n]). The remaining three objects appear in 
the H II region area. All of them are associated with a low value of the excitation indicator 
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(R < 0.5). Only one is located in the overlapping area between the PNe and the H II 
regions (MCMPc35) but it is indeed a low excitation object and thus presumably it is not 
a planetary nebula. Most objects without continuum were found to occupy the PNe region 
although they appear slightly displaced towards the upper border of the PNe band. Note, 
on the other hand, that the three objects where chemical abundances were determined 
(Section 4) lie inside the PNe area, and that all candidates in the PNe area show values 
of R from 0.3 to larger than 1.6, as expected for a normal population of PNe. In fact, 
only three candidates in the overlapping PNe and H n region area have 0.5 < R < 1 
(MCMP31, MCMP122, MCMP132). Although those objects can be misidentified with 
H II regions, they can equally well be low excitation PNe. The remaining candidates can 
be considered bona fide planetary nebulae. A total of ten candidates lie above the line 
plotted to mark the PNe region: seven in the H n region area and three in the SNR area. 
Three of the former, MCMP23, MCMP78 and MCMP127, have R < 0.5 and are probably 
unresolved H n regions. MCMP61 has higher excitation and its nature remains doubtful. 
The other three objects (MCMP18, MCMP25, MCMP28) have a high R and lie near 
the PNe locus where there are also numerous confirmed PNe. They are probably genuine 
PNe. MCMP41, MCMP42 and MCMP53 occupy the SNRs locus, each with different 
excitation. 

Concluding, no emission line object with some continuum emission (called 'unclas- 
sified ' in our previous work; MOO) can be considered to be PNe, whereas most of our 
candidate PNe (objects without continuum) are shown to occupy the same locus in the 
Ha/[N n] vs. Ha/[S n] diagnostic diagram as Galactic and extragalactic PNe. Three 
candidates which lie in the upper H n region area and have a low value of R are rejected 
as PNe, whereas four objects in that area showing high R are classified as PNe. Three 
candidates located in the SNRs region are rejected as PNe and another three candidates 
in the overlapping area between PNe and H n regions can be low excitation PNe or 
compact H n regions. All in all, 72% of our candidate PNe with good spectral data (i.e. 
with Ha, [N n], [S II]. [O III] lines well measured) are bona fide PNe; the remaining 28% 
will include H II regions, SNR and low excitation PNe. From the galaxy distribution of 
R discussed in our previous work (MOO), we did expect a contamination by H n regions 
of approximately 11%, a figure which is in rough agreement with the values found above. 
The probable nature of each object, as deduced from the preceding discussion, is specified 
in Table 1. 

3.3. Electron density and temperature 

The essential line for electron temperature determination, 4363A [O in], was measured 
in 6 PNe ( MCMP60, 65, 91, 93, 101, 122) . For the other 36 objects, it was not observed: 
for 21 objects it was too faint to be detected while in the remaining 15 the spectral range 
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did not include it (Section 2). The spectra of PNe MCMP91, 93, 101 and 122 had good 
signal to noise ratio, whereas MCMP60 and 65 were quite noisy. To test the quality 
of the spectra of the PNe where physical and chemical quantities can be computed, we 
compared the observed and theoretical values for the line ratios set by atomic constants 
and present in our spectral range: [O in] I(5007)/I(4959) and [N n] I(6584)/I(6548). The 
observed mean for the [O ill] ratio is 2.85± 0.40 and for the [N n] ratio is 2.85± 0.50, 
considering the sample of six nebulae, and 3.01±0.10 and 3.02±0.30 for MCMP91, 93, 
101 and 122, whereas the theoretical value, corrected for the energy of the levels, are 
respectively 2.96 and 2.98 (Osterbrock 1989f ■ We computed electron temperature and 
density for the whole six nebulae sample, avoiding chemical abundance determinations 
for the two PNe whose spectra were noisier, as deduced from their observed [O in] and 
[N n] line ratios and from their overall low S/N and taking out also the object MCMP122 
whose classification as a PN was doubtful. 

Electron temperatures and densities have been derived using the ratios [O ill] 
1(5007) /I(4363) for T e and [S n] I(6717)/I(6731) for n e . The most relevant source of 
error for the temperature-sensitive ratio is the error associated with the weak [O in] 
4363A line, whereas the error on the density ratio is based on the accuracy of the fluxes 
of both lines of the [S n] doublet. The uncertainty on T e is approximately 15% for the 
good S/N spectra, and somewhat higher for the remaining ones. Errors in the electron 
density are larger but density essentially does not enter into the chemical abundance 
determinations. We adopted T e [N li]=T e [0 in], i.e. the same temperature for the high 
and low ionization ions, because the [N il] 5755A line was not detected in our spectra 
and consequently T e [N II] could not be determined directly from the measurement of the 
[N n] I(5755)/I(6583) ratio. The nebular electron temperatures and densities for the 6 
PNe where we could measure the [O in] 4363A line are presented in TableEl In this Table 
their distances in kpc from the centre of M 33, their magjoni] (from MOO, extinction cor- 
rected with the c^s of Table and their de-reddened excitation parameter R— g^r^nj 
are also indicated. 



3.4. Chemical abundances 

As mentioned before, chemical abundances could be derived with a reasonable accuracy 
only for three out of the six PNe, namely MCPM91, 93, and 101, in whose spectra the 
[O m] A4363A was detected since their S/N was good also for the weakest lines. 

Helium ionic abundance have been computed considering the Case B recombination 
and using the effective recombination coefficients from Hummer & Storey l|1987fl for 
H/3 and He 2+ and from Brocklehurst l|1971fl for He + . The He I line strengths were cor- 
rected for the effect of collisional population of their upper state following Clegg 1|1987[) . 
The contribution of this effect is negligible in the case of low densities. The ionic metal 
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Table 2. Electron densities and temperatures. 





MCMP60 


MCMP 65 


MCMP 91 


n e (cm 3 ) 


300 


700 


1900 


T[onr] (K) 


12600 


14900 


11100 


d (kpc) 


3.1 


1.8 


4.3 


mag [0 iii] 


22.50 


21.98 


20.78 


R 


3.2 


3.2 


3.2 




MCMP93 


MCMP101 


MCMP122 


n e (cm -3 ) 


850 


100 


400 


T[onr] (K) 


9900 


11500 


13200 


d (kpc) 


3.8 


3.4 


0.8 


mag[oni] 


21.64 


21.19 


21.30 


R 


3.2 


4.1 


0.5 



abundances were derived from their collisionally excited lines solving the equations of 
statistical equilibrium. In order to calculate the total abundances, the unseen stages 
of ionization were accounted for using the standard ionization correction factors (ICF) 
(KB94). We used the same atomic data as KB94. 

We estimated errors in the chemical abundances of these three objects on the bases of the 
error propagation across the procedure of evaluation of abundances, taking into account 
the uncertainties in the observed fluxes of the relevant lines and in the electron tem- 
peratures. Errors due to the uncertainties in the atomic quantities were not considered 
due to the extreme complexity of that, also regarding obvious compensation of these 
errors among the various atomic quantities. We obtained formal errors of: 20-30% for 
He, 25-50% for O, 40-70% for N, and 30-70% for Ar. Errors are not given for S since 
its abundance is uncertain and it will be discussed in Section 4.2. The derived chemical 
abundances of the three PNe are listed in Table 

4. Discussion on the chemical abundances 

4.1. Comparison with other galaxies 

The median abundances of Helium, Nitrogen, oxygen, sulphur and argon of the three PNe 
are given in Table^J He abundance is given as He/H, whereas (N, O, S, Ar) abundances 
are in the usual log[^] + 12 units. In Table0]we show for comparison chemical abundances 
of 15 PNe in M 31 (JC99), the Milky Way (from the sample of KB94), the LMC and the 
SMC (CleggHUSI), the Sagittarius dwarf galaxy (Walsh et al. I1997j> and the elliptical 
galaxy Centaurus A (NGC 5128, Walsh et al. I1999[> , These data include all the most 
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Table 3. Chemical abundances of three PNe in M 33. 
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Ion/Elem. 


MCMP91 


MCMP93 


MCMP101 


He+/H 


0.17 


0.09 


0.14 


He 2+ /H 


0.02 


0.01 


- 


He/H 


0.19 


0.10 


0.14 


O + /Hxl0 4 


1.44 


2.37 


2.96 


O 2+ /Hxl0 4 


3.25 


4.20 


3.03 


icf 


1.08 


1.05 


1.00 


O/HxlO 4 


5.06 


6.92 


5.72 


N+/Hxl0 5 


1.90 


1.22 


0.77 


icf 


3.51 


2.92 


2.12 


N/HxlO 4 


0.67 


0.30 


0.16 


Ar 2+ /Hxl0 6 


1.49 


1.51 


1.48 


icf 


1.40 


1.52 


1.89 


Ar/HxlO 6 


2.08 


2.29 


2.80 


S+/Hxl0 6 


0.49 


0.85 


0.36 


icf 


1.16 


1.12 


1.06 


S/HxlO 6 


0.57 


0.92 


0.38 



recent information concerning chemical abundances of extragalactic PNe (Ford et al 
2002J. We also present the abundances of several H II regions: of M 33 (Kwitter & Allcr 
IP35T1 Vflchez et al. H"U55j) . of M 31 (Blair et al. HM^ . of LMC and SMC fDufour liUgljl 
and of the Orion nebula (Rubin et al. I1991f) . Finally, the figures can be compared with 
the solar abundances (Clegg©92). 

The median helium abundance of our PNe is close to that of the PNe in M 31 (JC99), 
and slightly higher than the Milky Way PNe one. The oxygen and argon abundances are 
in agreement with the values obtained by KB94. They are typically lower than the solar 
values, as observed in Galactic PNe. The nitrogen and sulphur abundances in M33 PNe 
are unexpectedly low. 

A first comment comes from the comparison of PNe and H II region chemical abun- 
dances (see Tabled}. The stellar evolution theory predicts that PNe should have enhanced 
abundances of He, N and C whereas the elements heavier than N should have abundances 
not appreciably different from those at the moment of the formation of their progenitor 
stars (cf. Iben & Renzini 1983). 

The comparison of abundances in PNe and H II regions for the galaxies of Table 4 
shows that in our Galaxy, LMC and SMC the theoretical predictions in all elements are 
verified. The prediction appears to be verified for He and N in M33 and for N in M31 
(note that in M31 the H n region abundance of He is not reported) whereas for O, S, 
and Ar the numbers do not support the predictions. Since the abundances both of PNe 



14 Magrini et al.: Spectroscopy of planetary nebulae in M 33 

Table 4. M 33 PNc chemical abundances in comparison with other abundances, with 
(N, 0, S, Ar) expressed in the usual lg[j|] + 12 units. 





Galaxies 


He/H 


N/H 


O/H 


S/H 


Ar/H 




M 33 (a) 


0.14 


7.50 


8.76 


5.76 


6.36 




M 31 (6) 


0.13 


8.03 


8.35 


7.38 


6.00 




Galaxy (c) 


0.115 


8.35 


8.68 


6.92 


6.39 


PNe 


LMC (d) 


0.105 


8.07 


8.44 


6.51 


5.93 




SMC (d) 


0.107 


7.84 


8.24 


6.48 


5.62 




Sgr (e) 


0.106 


7.41 


8.30 


6.40 


5.82 




Cen A (/) 


- 


8.02 


8.39 


7.10 


- 




M 33 (s) 


0.097 


6.98 


8.13 


6.81 


5.85 




M 33 (h) 


0.083 


7.34 


8.55 


6.95 


- 




M 31 w 




7.63 


8.86 


7.11 




H II 


Orion (j) 


0.100 


7.83 


8.60 


6.93 


6.65 


regs. 


Galaxy (fc) 


0.100 


7.57 


8.70 


7.06 


6.42 




LMC (fc) 


0.085 


6.97 


8.43 


6.85 


6.20 




SMC (fc) 


0.080 


6.46 


8.02 


6.49 


5.78 




Solar (d) 


0.098 


8.00 


8.93 


7.21 


6.56 



(a) This work (median values of 3 PNe) (t) Jacoby & Ciardullo ©99.3, (c) Kingsburgh & Barlow 
fT59H . (d) Clegg {15551 (e) Walsh et al. fiWfy . (/) Walsh et al. fTMfl l. (s) Kwitter et al. flMtl . 
(h) Vflchez et al fTM& , w Blair et al. JT5H51 . (j) Rubin et al. iTHMll . (fc) Dufour f[Sgfl . 



and H II regions are better measured in our Galaxy and the MCs than in M31 and M33, 
we infer from the above that the abundances of O, S, Ar in M31 and M33 might not 
be accurate enough to claim agreement or disagreement with the theory as far as the 
behaviour in PNe versus H n regions is concerned. We also compare the N/O of the three 
PNe in M33 to that found in H II regions in the same galaxy. According to Table 0] N/O 
is, on average, 0.05 for the three PNe, and similar to the N/O ratios of 0.07 and 0.06 
found by Kwitter et al. 119811 and Vilchez et al. 119881 for H II regions in M33. For any 
other galaxy in Tabled] N/O for PNe is much larger than for H II regions, showing the N 
enrichment occurring in the atmospheres of intermediate mass stars. It seems that such 
an enrichment has not occurred for the three PNe in M33. 

In the following, we first consider the abundance of S, which is unexpectedly low both 
in comparison with that in the H II regions of M33 and relative to PNc in our Galaxy. 
Next we address the abundance of N, which is higher than the corresponding abundance 
in the H n regions of M33, but lower than PNe in the Galaxy. 
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4.2. The sulphur abundance 

The unexpected low sulphur abundance in the PNe observed in M 33 can be explained 
by the fact that we are not observing the [S in] lines. In high excitation objects, such as 
PNe, most of the sulphur is twice ionized, and abundances derived using only [S n] lines 
are not representative of its total abundance, even if corrected with the nominal ICF. 
Thus we must consider our value as a lower limit to the true abundance of S. 

4.3. The nitrogen abundance 

The nitrogen abundance that we found is quite low for a galaxy of normal metallicity, as 
M 33 is, if the abundances of the (best observed) elements like He and O are compared 
with the corresponding values for the Galactic PNe. Our N/H median is 7.5, which is 
0.85 dex below the Galactic value, while the O/H abundances are comparable (see Table 
0J). We have investigated various possible reasons of this "under-abundance" . First, our 
chemical abundance determinations have been done using a single electron temperature 
(T e [0 III]) for all the ions. This approximation introduces errors in the calculation, es- 
pecially for the abundances of low ionization ions and, in turn, of the ICFs. oxygen and 
argon are dominated by high ionization states, while, in the optical, nitrogen can be ob- 
served only in its N + state; thus the use of T e [N II]=T e [0 III] affects the determination 
of N abundances. We have tried to test this effect assuming two different electron temper- 
atures. The low ionization temperature is usually derived from the [N n] I(5755)/I(6583) 
ratio, which provides T e [N n]. As the [N n] 5755 A line is not observed in our spectra 
of M 33, we tried to estimate whether a relation between temperature from [O in] and 
from the [N n] line ratio could be established in well studied Galactic PNe. For that 
we have considered the extensive work of KB94. We have found that, on average, an 
approximately linear relation holds between the [O in] and [N n] temperatures. Using an 
error-weighted least-squares fit, we found: 

T e [N II] = 0.23T e [O III] + 0.71 (1) 

where T e [0 in] and T e [N n] are expressed in 10 4 K. With this relation T e [N n] can be 
crudely estimated for the PNe with known T e [O in] . We computed chemical abundances 
both with T e [N n]=T e [0 ill] and with T e [N n] from Q. With T e [N II] < T e [0 III], the 
ionic abundance of N + is enhanced while the ICF(N) is lowered. The total abundance 
of N, however, does not change considerably so that the appreciable nitrogen deficiency 
remains. 

Another effect that we considered in the analysis of the possible N/H deficiency is the 
accuracy of the flux measurements of the [O n] doublet at 7325A. The most critical zone 
of our spectra is indeed the 7200-7400 A region, where the sky subtraction is particularly 
difficult. The [O n] flux is then quite uncertain. This implies possible errors on the + 
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abundance, in the ICF(N) and consequently on the N abundance. The oxygen abundance 
is dominated by the [O in] flux and the errors on + ionic abundance do not affect it 
greatly, whereas the ICF(N) depends on the ratio between the O and + abundances 
(see the prescriptions by KB94). We first consider the extreme case in which the sky 
subtraction had been plainly wrong in all three PNe and our mean ionic abundance 
of + was only one half of the value that follows from Table 3. The N abundances 
would turn to be comparable with the mean Galactic PNe N/H value. However, if the 
error in the sky subtraction was not so large and systematic the possibility remains of a 
real under- abundance in N for PNe of M33. To better evaluate this possibility we have 
compared the behaviour of versus for our PNe, which are disc PNe, with that 
of the Galactic sample of non-type 1 PNe of KB94. Our PNe fall in the lower-left part 
of the diagram, corresponding to the lower values of the PNe of the KB94 sample. 
Our |f , and also ^ are thus consistent with the corresponding values of the lower 
PNe of the KB94 Galactic sample. If our three PNe would be representative of the 
M33 PNe, their low ^ would indicate a poor N enrichment in the atmospheres of the 
progenitors of the PNe, and consequently a N under-abundance. 

Obviously, this under-abundance cannot be argued from data of only 3 PNe where 
the T e [N II] is not directly measured, but at the same time it cannot be excluded with 
the present data. On the other hand, it seems there is not any evident under-abundance 
of N in H II regions of M33. Their N/O ratio is equal, within uncertainties, to that of all 
the galaxies in Table 0] (except the determinations by Rubin et al. 119911 for Orion). The 
apparently low £ ratio in the studied PNe of M33 clearly deserves further investigation. 



5. Conclusions 

We have presented optical spectra of a sample of candidate PNe in the nearby spiral 
galaxy M 33. The location of the objects in the Ha/[S n] vs Ha/[N n] diagnostic diagram 
and their excitation parameter R = j^fMjj] allowed us to confirm that 26 of the 
observed 36 candidate PNe with reliable spectra are very likely genuine PNe, while seven 
other objects are likely compact H n regions and another three might be SNR candidates. 
None of the emission objects that were classified as "uncertain" by MOO (because having 
a non-negligible continuum emission) seem to be PNe: three are likely H n regions and 
two SNRs. Thus, these spectra confirm the overall reliability of the criteria adopted by 
Magrini et al. 120001 1200TEI 12002^1 to select candidate PNe in galaxies of the Local Group 
from ground-based narrow-band imaging. The degree of contamination found for M33 
PNe is approximately 20% and 10% for unresolved H n regions and SNRs, respectively. 

This work also confirms the potentiality of the use of PNe for the determination of 
chemical abundances for the old and intermediate populations of galaxies, as an alterna- 
tive to integrated light spectroscopy. With a 4 m-class telescope one can reach the limits 
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Fig. 4. vs N + /H: empty dots represent non-type 1 PNe from KB94 and filled dots 
the 3 PNe in M33. 
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of the Local Group and obtain a direct abundance determination for PNe, observing the 
relevant emission lines for density, temperature and metallicity analysis. We present a 
detailed physico-chemical analysis of the three best observed PNe in M33, finding that 
Hc/H, O/H and Ar/H abundances are in agreement with the Milky Way PNe values 
(KB94). N/H is unexpectedly low. This might be due to the weakness of the [O n] lines 
which are strongly involved in the calculation of the correction factors needed to deter- 
mine the total Nitrogen abundance from single low ionization ions. On the other hand, 
the Galactic behaviour of vs. ^ indicates that the value of ^J- for the PNe in M33 
is comparable with that of Galactic PNe with low nitrogen abundance. We cannot then 
fully exclude the possibility of a real nitrogen deficiency in PNe of M 33. 
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